KCL-PH-TH/2013-09 



Galaxy Rotation Curves in Covariant Hofava-Lifshitz Gravity 

J. Alexandre, M. Kostacinska 

j ean . alexandre@kcl .ac.uk, marty na. kost acinska@kcl .ac.uk 
King's College London, Department of Physics, WC2R 2LS, UK 



m 
O 

co 



X 



Abstract 



Oh- Using the multiplicity of solutions for the projectable case of the covariant extension of 
Hofava-Lifshitz Gravity, we show that an appropriate choice for the auxiliary field allows 
for an effective description of galaxy rotation curves. This description is based on static and 
spherically symmetric solutions of covariant Hofava-Lifshitz Gravity and does not require 
> • Dark Matter. 
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1 Introduction 

Although General Relativity (GR) has been accurately checked in the solar system, the 
small distance behaviour of gravity is still not understood, and several alternatives to GR 
have been proposed. Amongst these modified theories is Hofava-Lifshitz (HL) gravity [1], 
which is based on the Lifshitz approach and consists of introducing an anisotropy between 
space and time: when space is rescaled as x — )■ bx, time is rescaled as t — )■ b z t, where z is 
an integer. In Lifshitz models, Lorentz invariance must be recovered for z — 1, and the 
situation 2 ^ 1 leads to new renormalizable interactions: the convergence of loop integrals 
is improved by inclusion of higher orders in space derivatives, without introducing ghosts, 
since the number of time derivatives remains the minimum. A review of Lifshitz-type 
quantum field theories in flat space time can be found in [2], and the HL power-counting 
renormalizable theory of gravity is reviewed in [3]. 

A fundamental problem of the original model of HL gravity is the existence of an additional 
scalar degree of freedom for the metric, which can be understood as a Goldstone mode 
arising from breaking of 4-dimensional diffeomorphism. A solution to this problem has been 
proposed in [4] , where an auxiliary field A is introduced, such that its "equation of motion" 
leads to an additional constraint, eliminating the unwanted scalar degree of freedom of the 
metric. The resulting theory is invariant under a new Abelian gauge symmetry Us(l) 
which involves the metric components, the auxiliary field A and an additional auxiliary 
field v. This gauge symmetry, together with the 3-dimensional diffeomorphism of HL 
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gravity, can be shown to be equivalent to a 4-dimensional diffeomorphism at the lowest 
order in a post-Newtonian expansion, showing the equivalence with GR at long distances. 
For this reason, this modified version of HL gravity is called covariant HL gravity. We note 
however, that the long-distance limit is not obviously recovered: it has been shown in [5] 
that the equivalence principle is not automatically retrieved in the infrared, and that the 
meaning of the auxiliary field A and its coupling to matter are still open questions. 
Because of the anisotropy between space and time, HL gravity is naturally described in 
terms of the Arnowitt-Deser-Misner (ADM) decomposition of the metric, which expresses a 
space-time foliation. An important consequences of space-time anisotropy is the possibility 
of imposing the lapse function N to depend on time only. This situation is called the 
projectable case and leads to an interesting feature: the solutions of the equations of 
motion are not unique. Indeed, the Hamiltonian constraint, obtained by variation of the 
action with respect to N, leads to an integral equation, which does not have a unique 
solution, as will be seen in the present article. This multiplicity of solutions is independent 
of the above mentioned new gauge symmetry U%(1), and the different solutions for A in the 
projectable case belong to different gauge orbits. On the other hand, in the non-projectable 
case where N depends on both space and time, the Hamiltonian constraint leads to a 
differential equation, which has a unique solution, after fixing the constants of integration. 
The non-projectable case has been studied in [6] for static spherically symmetric solutions 
of covariant HL gravity. 

In this article we use the freedom of choice of the auxiliary field A in the projectable case to 
study the possibility of fitting galaxy rotations curves without the need for Dark Matter. 
One can argue that the ambiguity in the choice of A leaves its physical interpretation 
unclear. This may be understood as one of the incomplete aspects of covariant HL gravity, 
and one could conjecture the existence of an additional symmetry which would restrict the 
set of solutions for the auxiliary field A, but this still remains to be studied. We note that 
the multiplicity of solutions, arising from a Hamiltonian constraint expressed in terms of 
an integral over space, has been discussed in [7] as an alternative to Dark Matter, in the 
context of the original HL gravity. Fits to galaxy rotation curves have been studied in the 
framework of the original HL gravity [8], where a detailed comparison with experimental 
data is presented. Also, an interesting analogy between HL gravity for z = and Modified 
Newtonian Dynamics theories has been described in [9]. 

In section 2 we review the static and spherically symmetric solutions of covariant HL 
gravity, and in section 3 we show how galaxy rotation curves can be described by the 
covariant extension of HL gravity. Although the curves we use were derived from Dark 
Matter models, we use them as experimental data fits. Finally, we conclude by discussing 
the possiblity to describe the solar system with our solution. 
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2 Static and spherically symmetric solutions of co- 
variant HL gravity 

The static and spherically symmetric solutions of covariant HL gravity were derived in [13], 
and we review here the aspects relevant to our present study. 

2.1 Action 

The ADM metric we consider is 

ds 2 = -c 2 N 2 dt 2 + 9ij (dx l + N*dt) (dx j + N j dt) , (1) 

where N and N l are the lapse and shift functions respectively, and ^ is the three- 
dimensional space metric. We consider the anisotropic scaling with z — 3, for which 
the different operators, allowed for the model to be power-counting renormalizable, have 
the mass dimension of 6, at most. We note that in this case the mass dimension of the 
speed of light is [c] = z — 2 = 2, and we will keep c in the different expression for the sake 
of clarity. 

The covariant Horava-Lifshitz action, in the absence of cosmo logical constant, is 

S = J dtd 3 x^ (iV [KijK ij -K 2 -V + u& ij (2K ij + ViV»] - AR (3) ) , (2) 

where g is the determinant of the three-dimensional metric gij, with curvature scalar R( 3 \ 
A is an auxiliary field (without kinetic term) of mass dimension z + 1 = 4 and v is an 
auxiliary field of mass dimension z — 2 = 1. In the above expression, the extrinsic curvature 
is 

K H = ^ i9ij - ViNj - VjNA , i, j = 1, 2, 3 , (3) 
where a dot denotes a time derivative, and 

Qij = ^ij _ \ R $)gij ^ (4) 

2 

where where R^ lj is the Ricci tensor corresponding to the three-dimensional space metric 
g^. Note that a more general kinetic term K^K^ — XK 2 can be considered, with A a generic 
parameter. The latter should be equal to 1 if one wishes to recover GR in the infrared, but 
in the framework of HL gravity, no symmetry imposes this parameter to be equal to 1. It has 
been shown that covariant HL gravity with A ^ 1 leads to inconsistencies, if one compares 
the predictions of the model with solar system tests [10], and we therefore focus here on the 
situation where A = 1 only. Independently of these tests, the mathematical consistency of 
covariant HL gravity for A ^ 1 has been studied in [11], and the corresponding Hamiltonian 
structure in [12]. 

The potential term V contains up to six spatial derivatives of the metric g^: 

V = _ c 2^(3) - ai {R^) 2 -a 2 R^Rf -^{R^f-foR^R^Rf (5) 
-^Rf )l Rf )3 Rf )k - (3 4 R^V 2 R^ - hVRffriR®* , (6) 
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where the mass dimensions of the different couplings are 

M = H =z-l = 2 , [ft] = l(3 2 ] = [f3 3 ] = [ft] = [ft] = z - 3 = . (7) 
We consider the most general, static, spherically symmetric metric, of the form: 

ds 2 = -c 2 N 2 dt 2 + j^(dr + n(r)dtf + r 2 (sin 2 6d6 2 + d<f) 2 ), (8) 

where n(r) = N r (r) is the radial component of shift function, and N r = n(r)/f(r) since 
g rr = l//(r). Note that n has mass dimension z — 1 = 2. 

2.2 Constraints and equations of motion 

The variation of the action (2) with respect to the different degrees of freedom leads to 
the following constraints, or equations of motion, where a prime denotes a derivative with 
respect to the radial coordinate r. 

• The variation with respect to A gives = 0, or equivalently 

r/' + /-l = 0, (9) 
which imposes the function f(r) in the metric (8) to be 

f(r) = 1 - ™ , (10) 
where B is a constant of integration. 

• The variation with respect to v gives 

e ij ViVjU + Q^Kij = 0, (11) 

In what follows we will assume the Gauge fixing of v — 0, then the above constraint 
gives Q lj Kij = 0, which is satisfied for spherically symmetric solutions. 

• The variation with respect to iV gives the Hamiltonian constraint 

poo -r>2 

/ dr^==(K lJ K^-K 2 + V) = 0, (12) 

which is an space-integral equation since N depends on time only. 

• The variation with respect to n gives 

f(r)n(r) = , (13) 
such that either / is a constant (B = 0) or n = 0. 

• The variation with respect to / gives 

where the time has been rescaled such that N — 1. 
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2.3 Solutions 



The different solutions of the previous set of equations were derived in [13] and here we 
shortly review the different cases. 

• n = : the solution of the differential equation (14) gives short-distance corrections 
to the Schwarzschild solution, details of which can be found in [13]; 

• n^O and A = 0: leads to the Schwarzschild solution in the Painleve-Gullstrand co- 
ordinates (n oc r~ 1//2 ), and the Hamiltonian constraint is automatically satisfied; 

• n^O and A ^ 0: corresponds to the situation with multiple solutions for A and n, 
on which this article focuses. Since B = 0, we have / = 1 and V = 0, such that the 
equations (12) and (14) give respectively 

/ drirn 2 )' = 
Jo 

rA' + 2(ra 2 )' = , (15) 



with solutions 

r,' 2 (r) - 

r 2 w 2r Jo 
where C is a constant of integration and the auxiliary field A satisfies 



>r(r) - \a{t) • 1 / dp A{p) , (16) 



L 



drrA'ir) = . (17) 

o 

The flexibility in the choice of solution for A will help describe galaxy rotation curves, 
as explained in the next section. 



3 Galaxy rotation curves 

In what follows, we consider each star in the spiral arms of the galaxy as a test particle, 
moving on a circular trajectory under the influence of the potential 4>(r), generated by the 
centre of the galaxy, which is assumed to be static and spherically symmetric. Our approach 
is to start from a stellar velocity distribution and derive the corresponding expression for 
the auxiliary field A, in the above situation with n^0 and A ^ , where multiple solutions 
for A are allowed. We note the following few points: 

(i) We consider vacuum solutions of covariant HL gravity, for which one cannot in principle 
describe the region in the centre of the galaxy. Nevertheless, as we will see in case of the 
following velocity profiles, a consistent solution for A can be found for < r < R, where 
R is a typical radial length describing the galaxy. To be consistent with the usual studies 
of galaxy rotations, we will choose here R = r 2 oo, corresponding to the virial radius of 
the galaxy. The latter is defined as the distance from the centre of the galaxy, where the 
density p 2 oo is 200 times the critical density p c of the Universe 

p 200 = 200 Pc = 200-— , (18) 
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where H is the Hubble constant; 

(ii) For r > r 20 o, we assume the solution A = 0, which leads to the usual Schwarzschild 
solution outside the galaxy. Indeed, far from the galaxy, one expects to see the Newtonian 
potential. This solution implies that the Hamiltonian constraint is expressed in terms of 
an integral over the finite range [0, r2oo] of radial coordinate r, which we can impose to 
vanish by fixing the constants of integration; 

(Hi) For r < r 20 o, since the shift A — > A+constant does not have a physical implication, 
one can expect to find a consistent solution for A which is continuous at r = r 20 o, i- e - 
^4( r 2oo) = 0. On the other hand, one cannot expect the auxiliary field to be differentiable 
at r = r 2 oo, but this is not necessary for the consistency of covariant HL gravity: the 
second derivative of A does not appear in any equation of motion or constraint, so the first 
derivative can be discontinuous. 



3.1 Gravitational potential 

On a circular trajectory, the relation between the speed v of a star and the Newtonian 
potential is 

- = (W) 

r ar 

where is obtained from the g 00 component of the low energy effective theory 

_c 2 + n 2 = _ c 2 (i + 20) , (20) 

and is thus given by 

0(0 = • (21) 

We note here that, in the situation where the second auxiliary field v does not vanish, the 
Newtonian potential is instead given by 

m = _ ("W-v^(r))' (22) 

and is independent of the Us(l) gauge choice for v and A. Also, is dimensionless, so 
that the speed v we use in this article is also dimensionless, which corresponds to a usual 
definition of speed, for isotropic space-time. The corresponding "Lifshitz velocity" is cv, 
with mass dimension z — 1 = 2. 

From the relations (16) and (19), the speed v and the auxiliary field A are then related by 
rA'(r) d ( r r , v 2 (p)\ . r r , v 2 (p) 

= * ( r / d '^r j = - V) + / d '^r • < 23 > 

and we have to check that the Hamiltonian constraint (17) is satisfied, in order to show 
the consistency of the approach. 

We note from eq.(23) that, for r > r 20 o, the solution A = gives 

v 2 (r)+ [ r dp^ = 0, (24) 
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which leads to the Newtonian velocity v oc r 1 / 2 , and the non-trivial solutions for A which 
are derived below are valid only inside the galaxy disc, with r < r 2 oo- 



3.2 Navarro, Frenk and White profile 



This profile originates from Cold Dark Matter halo models and is characterised by the 
"cusp" shape of the density distribution. Its validity has been tested on variety of obser- 
vational results, including Low Surface Brightness galaxies in [14] and spherical galaxies 
and clusters in [15]. For our purposes, we use the circular velocity profile derived from the 
mass density profile in [16], based on the assumption of a Dark Matter halo and ignoring 
the contribution of baryons, 



v(r) = v J In 1 H 

V r V r 200 



/ 1 + ar/r 200 ' 



where a is the concentration parameter and 



V200 



v /ln(l + a)-a/(l + a) 



(25) 



(26) 



where t>2oo, the virial velocity, is the circular velocity at the virial radius. Some typical 
values are [17] 



Galaxy 


a 


v 200 [km/s] 


^200 [kpc] 


NGC 2403 


10.9 ±0.6 


106.1 ± 1.9 


88.4 ±2.6 


NGC 3198 


11.2±0.43 


104.0 ±0.7 


86.6 ± 1.6 


NGC 3521 


14.0±12.6 


122.5 ±20.4 


102.0 ± 18 



With the profile (25), eq.(23) leads to 

rA'(r) a 
4c 2 Vq l + ar/r 2 oo 

where 4>i is a constant of integration, and therefore 

A(r) 



(27) 



Ac 2 Vq 



(0i - a) In + a In ( 1 + a-^— ] ± A 1 , 

Vr 20 o/ V r 20 o/ 



(28) 



where Ai is a constant of integration. Since we consider the solution A(r) = for r > r 2 oo 
and impose A(r) to be continuous at r = r2oo, then A\ = — aln(l ± a) and 



Mr) 
A(r) 



(0i — a) In 

V r 200 

if r > r 20 o 



± a In 



1 + ar/r-200 
1 + a 



if < r < r 20 o 



(29) 
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Finally, the Hamiltonian constraint 

fOO /T200 / (7 \ 

° = /„ *rrm =*c\lj o dr ^ 1 - rT - 7 _j , (30) 

is satisfied if 



0i = ln(l + a) . (31) 

We note that the solution (29) is consistent for any relevant values of the parameters 
a, f ,r 2 oo and therefore allows the description of a whole range of galaxies. 



3.3 Pseudoisothermal profile 

The pseudoisothermal mass density profile assumes an existence of a "cored" dark matter 
halo component in the galaxy [18] , with mass density of an approximately constant value 
in the central region of the galaxy, for r < R c . Empirically motivated, the model is often 
contrasted with the aforementioned NFW "cuspy" profile, while evidence shows that it 
provides better fit of the galaxy rotation velocities [14]. The velocity profile is 



v(r) = v Jl — — arctan ) , (32) 



r \R 



where v is given by 

v = ^4irGp R 2 c , (33) 
and po, the central density, is the density within Rc- Typical values are [17] 



Galaxy 


R c [kpc] 


po [lO- 3 M pc- 3 ] 


^200 [kpc] 


NGC 2403 


2.51 ±0.32 


59.1 ± 14.3 


88.4 ±2.6 


NGC 2841 


1.36 ±0.75 


674.8 ± 736.4 


159.0 ±3.7 


NGC 3621 


5.88 ±0.32 


13.0 ± 1.1 


106.6 ±4.0 



We note that in the case of NCG 2403, both the NFW and the present profile can be shown 
to produce fits of comparably good quality. 
With the profile (32), eq.(23) leads to 

where 2 is a constant of integration. The same steps as those describe for the previous 
profile lead to the solution 

^4 = 2 lnf— ) + i [ ,Rcd i\n{l + t 2 ) + A 2 if 0<r<r 200 (35) 
Ac 2 v^ \r 20 o/ 2 Jo t 

A(r) = if r > r 200 , 
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where the constants of integration are 

A 2 = --J Q T ln(l + t 2 ), (36) 
for the continuity of A(r) at r = r 20 o, and 

1,.. , r 200 \ , 1 Rc /' '''21111 



for the Hamiltonian constraint to be satisfied. 



4 Conclusion 

The multiplicity of solutions in covariant HL gravity has been used to describe galaxy 
rotation curves, without the need for Dark Matter. This multiplicity is a consequence 
of the breaking of 4-dimensional diffeomorphism, in the projectable case, where the lapse 
function is imposed to depend on time only. We considered two specific velocity profiles 
and showed in each case that a corresponding consistent solution can be derived for the 
auxiliary field A. Clearly, similar results can be expected with other profiles. A more 
detailed analysis would consist in taking into account the finite star density at the centre 
of the galaxy, as well as its rotation. 

It is interesting to note that the solutions (29) and (35) actually do not only describe 
galaxy rotations curves, but also the usual Schwarzschild solutions for the Solar System, 
by taking the limit r2oo — > such that A = for all r. In this sense these solutions provide 
an effective unified description of gravitational effects, from the solar system scale to the 
galaxy scale, where local Lorentz violation appears gradually, as the observational length 
increases from the the size of the solar system to the size of the galaxy. 
We finally note that, although HL gravity was initially developed to modify ultraviolet 
behaviour of gravity, our present work uses another consequence of local Lorentz symmetry 
violation, which is the existence of a non-unique solution to the Hamiltonian constraint. 
It would be interesting to see if non-unique time- dependent solutions could also allow for 
an effective description of Cosmology, including an alternative to Dark Energy. 
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